Enhanced nonperturbative effects in jet distributions 
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We consider the triple differential distribution dT /dEjdm 2 jdQ,j for two-jet events at center of 
mass energy M, smeared over the endpoint region m 2 j M 2 , \2Ej — M\ ~ A, Aqcd < A < M. 
The leading nonperturbative correction, suppressed by Aqcd/ A, is given by the matrix element of 
a single operator. A similar analysis is performed for three jet events, and the generalization to any 
number of jets is discussed. At order Aqcd/A, non-perturbative effects in four or more jet events 
are completely determined in terms of two matrix elements which can be measured in two and three 
jet events. 



This paper studies the semi-inclusive two-jet distribu- 
tion AT / (dEjAm 2 jc\Vlj) in high energy processes such as 
Z — > qq, where Ej and mj are the energy and invari- 
ant mass of one of the jets, J. No restriction is placed 
on the kinematics of the second jet, J'. We are inter- 
ested in non-perturbative effects that are enhanced in 
the endpoint region \2Ej - M\ < M, m 2 , < M 2 . In the 
endpoint region the jet J is narrow, and has energy close 
to M/2. At leading order in the parton model descrip- 
tion, the quark and antiquark each have energy M/2, and 
irij = 0. Gluon radiation and non-perturbative effects 
such as hadronization can cause Ej and m 2 to deviate 
from their leading order parton model values. The aim 
of this letter is to characterize these non-perturbative ef- 
fects in QCD in terms of operator matrix elements and 
study the relationship between nonperturbative effects 
in two jets, and events with three or more jets. These 
non-perturbative effects do not depend on the precise jet 
definition; for example, a suitable choice is the Sterman- 
Weinberg prescription |]|. 

In the endpoint region, the decay distribution for two 
jet events can be written as 



dr 



dEjdm 2 jdn.j 



T(M,n,j)J(M, m 2 j/M) 



(1) 



x J dk+J(M, k+)S (2Ej - M - k + - m?j/ {2Ej)) . 

Here the hard function T and the jet function J are cal- 
culable perturbatively and S is a nonperturbative shape 
function. The precise form of the perturbative expan- 
sions of these functions depends on the jet definition 
used. In the endpoint region this distribution must be 
smeared over a suitable region of Ej and m 2 to be phys- 
ically meaningful. This general form for the decay distri- 
bution has been obtained Q previously using standard 
"diagrammatic" factorization methods Q • We begin this 
paper by reviewing the derivation of Eq. Q using soft 
collinear effective field theory (SCETV 

The effective field theory SCET 0,11 is appropriate 
for the kinematic region of interest. SCET describes the 
interaction of collinear and ultrasoft (usoft) degrees of 
freedom with momenta scaling as p c — (n ■ p,h ■ p, p ) ~ 
M(A 2 , 1, A) and p us ~ M(A 2 , A 2 , A 2 ). The light-like vec- 



tors n and n satisfy n° = n = 1 and n = — n and the 
perpendicular components of any four vector V are de- 
fined by V? = V - (n- V)n^ /2 - {n ■ V)n^/2. For our 
analysis, A 2 ~ A/M -c 1. The effective theory provides 
a simple method for the factorization of hard, collinear 
and usoft degrees of freedom at the operator level. For ex- 
ample, the factorization of usoft degrees of freedom arises 
because they can be decoupled from the collinear degrees 
of freedom using a simple field redefinition. SCET gives 
field theoretical definitions of the various ingredients in 
Eq. p. 

The first step in the SCET derivation of Eq. (JTJ is 
matching the full theory current onto SCET. The cur- 
rent in the effective theory at leading order in A is 



r 



(2) 



where = gvl±+9A"f±"f5, 9v,A are the vector and axial 
couplings of the quarks to the Z boson, and C(V\P) is 
the matching coefficient which is one at tree level. The 
field denotes a collinear fermion in the n direction, 
and we have used the convention 



p 



-ip-x 



(3) 



where p is the label momentum which contains compo- 
nents of order 1 and order A. The order A 2 compo- 
nents are associated with the spacetime dependence of 
the fields. 

The label operators V, V pick out the order one mo- 
menta of the collinear fields and W n (x) denotes a Wilson 
line of collinear gluons along the path in the n direction. 
The Wilson lines W n ,n are required to ensure gauge in- 
variance of the current in the effective theory [5j. The 
Lagrangian of the effective theory does not contain any 
direct coupling of collinear particles moving in the two 
separate directions defined by n and n 7], however they 
can still interact with one another via the emission of 
usoft gluons. The coupling of collinear and usoft degrees 
of freedom in the Lagrangian can be eliminated via the 
BPS field redefinition |(| 



A„ 



yt/i y 

j /in ± n. 



(4) 
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where 



Y n (z) = exp 



ig / dsn ■ A us (ns + z) 



(5) 



denotes a path-ordered Wilson line of usoft gluons in the 
n direction from s = to s = oo, since we are dealing 
with final state collinear fields. [For annihilation, Y n is 
from s = — oo to s = 0.] It is well known from the 
"diagrammatic" approach to factorization that the usoft 
degrees of freedom couple to collinear degrees of freedom 
via a Wilson line . 

To derive the factorized form in Eq. , we start from 
the general expression 

dr 



d 4 r 2M ^ K ' 

JX 



S i (M-pj-p x )S i {r~p J ) 



x|e M (J,X\f\0)\< 



(0) 
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where the sum includes the phase space integrations over 
all the particles in the final state, j M denotes the current 
producing the qq pair, and is the polarization of the 
decaying particle. The final state hadrons have been di- 
vided into those in the quark jet J, with total momentum 
r and the remaining hadrons (including the antiquark jet) 
which form X . 

Using these definitions, the matrix element of the cur- 
rent in Eq. JBJ becomes 

(JX\r\0) = g C( q +, fc") {j\T{[CnW n ]l a }\0) (T%, 

{[YnYX"[WlU^b}\o) , (7) 

where a, [3 (a, b) denote spin (color) indices, and the sub- 
script [-]q denotes the total label momentum of the oper- 
ator inside. In Eq. (J7J), and in the subsequent equations, 
we will use the notation n ■ p = p + , n ■ p = p~~ for any 
four-momentum p. 

The operators inside the matrix elements in Eq. Q are 
time-ordered. The operator T{F„} is a Wilson line where 
the operator time-ordering agrees with the path-ordering 
of matrix multiplication, T{Y^.} = Y n . The operator 
T{Yj}} has operator time-ordering in the opposite order 
as the path-ordering of matrix multiplication. The two 
orderings can be made to coincide by taking the trans- 
pose of all the matrix indices, so that T{Y^ a b } is the Wil- 
son line in the 3 representation, i.e., T\Y^ a b } = Y n b a . 
Similar results hold for the anti-time ordered products. 

Using Eq. J7J, the differential distribution becomes 
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FIG. 1: Pictorial representation of the Wilson lines occurring 
in the shape function S(k) defined in Eq. Ijlljl . 




The two terms in square bracket can be simplified, fol- 
lowing Ref. For the second term we have 



E 



j|r{[f fi w fl ]Ho)}| 

S a c Mr~), 



(9) 



where the last equality defines Jf(r ). 

Combining the first term with the <5" in Eq. gives 



E 

X 



X\T{[Y n Y^][W^ n ] s }\0 
= 3 j dk+J~ s {k+)S(-s+ - k+) (Pj . 
The shape function S(k + ) is defined by 



0pa 

(10) 



S{k) 



da 
2^ 



-,iku I q 



Tr T 



{[Y n Y^}(nu)]T{[YnY^}(0)} 



Yl c a Yl b S(k - in ■ d)Y nb e Y n \ 



(1 



The Wilson loop, shown in Fig.U is regulated to preserve 
reparametrization invariancc 9], and is a function of n ■ 
nu. S(k) is normalized so that /_ dk S(k) = 1. 

The jet function in Eq. © depends on f, which has 
both + and _L components. One can align n with the 
jet-axis, i.e. choose rj_ = 0. With this choice, we write 



J(r+ r~) = J r+ , 0± (r-). 



(12) 



Combining Eqs. i|5)l- HI2|) and using d 4 r = 
{M/4)dEjdm 2 jd£l,j gives Eq. QJ with the hard function 
defined by 



T(M,n.j) 



^j|C(M,M)|%e"*Ti: [T^fS] ■ (13) 



(8) 



The simple form for the shape function in Eq. (jl 1|) 
arises because the decay distribution is totally inclusive 
on the unobserved usoft degrees of freedom and on the 
other jet J'. This would not be the case for the thrust 
distribution near T = 1, since the deviation of T from 
unity depends on the momentum of the usoft degrees of 
freedom. 

We can simplify the discussion by restricting ourselves 
to the leading order in perturbation theory, where 



C(M, M) = 1, J(M, k+) = 2n5{k + ). 



(14) 
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Perturbative corrections to this result can be included 
by calculating the two jet functions J and the Wilson 
coefficient T to higher orders in perturbation theory. For 
Sterman- Weinberg jets pj, the perturbative corrections 
contain logarithms of the cone angle 6 and the minimum 
energy cut E c . In the effective theory, logarithms of 5 
arise for precisely the same reason as in the full theory; 
one has to separate three collinear-field final states into 
two-jet and three-jet events. Logarithms of E c arise in 
the effective theory from the separation into soft and 
collinear fields. One can sum these logarithms using 
renormalization group methods in the effective theory. 
Integrating Eq. over m 2 j using Eq. (|14|l gives 



dr 



dr(°) 

2^^S(2Ej - M) + 



dVL 



J 



(15) 



where dl^ ) /dfij is the parton model differential decay 
rate and the ellipsis denotes subdominant perturbative 
and power corrections. To this order the shape of the 
jet energy spectrum is determined entirely by the non- 
perturbative shape function. Non-perturbative usoft ra- 
diation is not preferentially directed along the jet direc- 
tions and will reduce the energy associated with the jets. 
Thus, S{k) = for k > 0. 

To proceed further we consider observables in which we 
smear the jet energy distribution over a region \2Ej — 
M\ ~ A <C M, with A > Aqcd |H|. In this re- 
gion there are enhanced non-perturbative corrections 
suppressed only by powers of Aqcd /A. Note that for 
the total rate the non-perturbative corrections are much 
smaller since they are suppressed by powers of Aqcd/-^- 
To be more precise, consider a smooth window func- 
tion wa(Ej) which when integrated over Ej is normal- 
ized to unity and has support concentrated in the region 
\Ej - M/2\ < A. Then, 



d^) A = I dEjWA{Ej) 



dr 



dEjdrij ' 



(16) 



is such an observable. For it we can expand the shape 
function in a power series, given by 



S(k) = S(k) - S'(k) (0|Oi|0> + ^S"(k) (0|O 2 |0) 
where 



O m = 3 Tr 



We can rewrite Oi as 



Y£(in-D) m Y n 

1 f°° 
Oi= 0\ = -Tr / dsQ nn (ns) 

Jo 



.(17) 



(18) 



where 



Qmn 2 {ns) = Y ii (ns,0) t n£ra 2 'G /11 ,lfi(ns,0) 



(20) 



Here, G^ v is the gluon field-strength (with a factor of the 
strong coupling absorbed into it) and Y fi (ns,0) denotes 
a usoft Wilson line along the n direction from to hs. 



Note that the Fourier transform of S(k) is a function only 
of the combination (n • n) u, so the matrix element of O m 
must have the form 



(0|O m |0) = (n-n) r 



A q 



(21) 



where A^ is a number of order Aq CD . If the observed jet 
is the antiquark jet, then in ■ d in Eq. (|19|) is replaced by 
in ■ d, and the leading correction is given by the vacuum 
matrix element of the antiquark operator 



Oi = ±TrFl(in-DrY n ], 



(22) 



which is (n ■ n) A q n . By charge conjugation invariance, 
A q = A q 

The nonperturbative corrections to the differential de- 
cay rate in Eq. (|15(l are singular at Ej = M/2. This is 
also the case for the perturbative corrections. However 
for the smeared observable in Eq. I|16|) the expansion of 
the shape function gives, 



ar 



J J A 



dflr 



l + w' A {M/2)A\ 



(23) 



Since the matrix elements of the operators scale as A^ ~ 
Aq CD and the m'th derivative of the window function 
(evaluated at Ej = M/2) scales as 1/A m , the square 
brackets on the right hand side of Eq. (|23|l contain an 
expansion in powers of Aqcd/ A. 

A similar calculation also applies for events with more 
than two jets in the final state. For example, three jet 
events contain an energetic gluon radiated off one of the 
quarks at large angle. The subscripts 1,2,3 refer to the 
quark, gluon and antiquark jets, respectively. One of 
the jets is unobserved, and is summed over. Consider 
the case where the antiquark jet (jet 3) is unobserved, 
and the quark and gluon jet are observed. The light like 
vectors ni,2 are defined by n% = (1, ni), ni = (1,112) 
where 11^2 are unit vectors in the direction of the two 
observed jets. One can the construct the third vector 



^ = 1 



n 3 



\Eim +E 2 n 2 \ ' 



(24) 



using only information from the two observed jets. A 
similar argument holds if the quark jet, or the gluon jet, 
is the unobserved jet. 

Let i and j denote the two observed jets, and r denote 
the unobserved jet, where (i,j,r) is some permutation of 
(1,2,3). We find at tree level that 



dr, 



dr 



(0) 



(!9) d/;,d<>,d/.,d<>, dEidflidflj 



(25) 



xSij (n r ■ niEi + n r ■ njEj — M) , 

for the differential decay distribution for the observed jets 
i and j in terms of a non-perturbative shape function Sij . 

In the three-jet case, the current in SCET contains the 
gluon field A n2 . Factoring the usoft degrees of freedom 
results in the shape function 



4 



SvW = \ /J e ife "Tr ^^{[y^T^Y^ . (26) 



r 



Thus, a different shape function determines the usoft 
physics in three-jet events, and in general a new non- 
perturbative function is required for each additional jet. 
This is not surprising since the color structure is very dif- 
ferent for events with different numbers of jets. The time- 
ordered product T{Y,l a b Y nc d } is equal to (1/3)5^ + 
ty n AB (T A ) c h (T B ) a d where y„ is the adjoint Wilson 
line from to oo in the n direction. 

In the kinematic region \M — Ein r ■ m — Ejn r -rij\ ~ A 
an expansion of iSy(fe) analogous to the one in Eq. i|17|) 
can be performed, and the first correction to the smeared 
rate is determined by the operator 



01 



- ^TrYl(in r -D)Y ni + ±TrYi 3 (in r -D)Y r , 



1 



+-Tr^t (i nr . D )y n 



(27) 



The matrix element of this operator is given by 



(0|Oi 3) |0) 



{n 1 +n 3 )A q 1 +n r -n 2 A 9 1 . (28) 



where A\ is the same number that occurs in two-jet 
events, and is defined by Eq. (|21|) and Af is a new number 
of order Aqcd ■ 

In three jet events with an unobserved gluon jet, 

(3) 

n r = n 2 and the vacuum matrix element of 0\ is 
n 2 ' {ni + ^3) A\, and is completely determined by the 
two jet case. If the unobserved jet is the quark jet, 



then n r = 



and the vacuum matrix element of O 



(3) 



is (ni • n 2 ) A\ + (ni ■ n 3 ) A\, and if the unobserved jet 



is the antiquark jet, the vacuum matrix element of 0\ 



(3) 



is (713 • n 2 ) A\ + (m ■ 713) A\. If one cannot distinguish 
quark and gluon jets, then the three jet distributions 
are given by averaging over the unobserved jet being the 
quark, gluon or antiquark jet. 

One can repeat the above analysis for events with four 
or more jets. In the endpoint regions, the generalized 
shape functions can be expanded as in Eq. (|17|l . and 
the leading non-perturbative corrections are expressed in 
terms of A\ and A\ which characterize two and three- 
jet events. An important feature is that, even though 
the shape functions for events with different numbers of 
jets are very different, at order A/Aqcd no new non- 
perturbative quantities enter for four or more jets after 
smearing by an amount Af > A > Aqcd over the end- 
point region. The reason is that the only Wilson lines 
that can arise from hard quark anti-quark and gluon ra- 
diation are in the 3, 3 and 8 representations, and these 
already occur in the two and three jet cases. One can also 
predict the distributions for two jets plus a hard photon 
in terms of A\ measured in two jet events, since the ad- 
joint Wilson line from gluon emission does not enter. 

A\ and A\ can be determined by fitting to the ex- 
perimental data on two and three-jet events. They can 
also be computed numerically by light-cone lattice gauge 
theory methods [T(i| , since they only involve Wilson lines 
along light-like directions. 
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